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TECHNICAL NOTE NO. 1706 


STAB IT, T TY DERIVATIVES OP TEEN RECTANGULAR STINGS 
AT SUPERSONIC SPEEDS 

WING DIAGONALS AHEAD OP TIP MACH LINES 
By Sidney M. Harmon 


SUMMARY 


Theoretical results are obtained, "by means of the linearized, theory, 
for the surface— velocity-potential functions, surface— pressure distribu- 
tions, and. stability derivatives for various motions at supersonic speeds 
of thin flat rectangular wings without dihedral. The investigation 
includes steady and accelerated vertical and longitudinal motions and 
steady rolling, yawing, sideslipping, and pitching for Mach numbers and 
aspect ratios greater than those for which the Maeh line from the leading 
edge of the tip section intersects the trailing edge of the opposite tip 
section. The stability derivatives are derived with respect to principal 
body axes and then transformed to a Bystem of stability axes. In the 
case of yawing, a treatment for the infinitely long wing which takes 
account of the spanwise variation in the stream Mach number is extended 
to the finite' wing, and a plausible, although not rigorous, solution is 
obtained for the wing tip effects. 

The results for this investigation showed that positive yawing at 
supersonic speeds may produce a negative rolling moment in contrast to 
the behavior at subsonic speeds where a positive rolling moment would be 
produced. The attainment of supersonic speed should produce a signifi- 
cant c hang e in the positive direction of the yawing moment per unit 
rolling velocity. The results also indicate that unstable tendencies 
are produced by vertical accelerations if 

a/m 2 - 1 > - — — 

3 

where A is wing aspect ratio and M is stream Mach number. 


INTRODUCTION 


Recent developments in supersonic airfoil theory (references 1 to 4) 
have led to the calculation of many of the supersonic stability derivatives 
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for various plan forms. In references 5 to 8, various theoretical super- 
sonic stability derivatives for small disturbances are presented for 
thin flat brings of delta plan form. In reference 9j the supersonic 
damping due to ro lling is given for triangular, trapezoidal, and related 
plan forms. 

In the present paper the methods of references h, 10, and 11, which 
are based on the linearized theory for a uniform stream Mach number, are 
used to derive the supersonic surface— velocity— potential functions for 
thin flat rectangular wings without dihedral in steady and, accelerated 
vertical motions and steady rolling, sideslipping, and pitching motions. 
The potential functions that are obtained are then used to derive formulas 
for the pressure distributions and the stability derivatives for the 
foregoing motions and also for steady yawing. In the case of yawing, a 
simple treatment given in reference 7 for the infinitely long wing, which 
takes account of the spanwise variation in stream Mach number associated 
with yawing, has shown that the assumption of a uniform Mach number is 
far from adequate to describe the compressibility effects. This treat- 
ment is extended herein in order to evaluate the wing tip effects for the 
yawing finite-span wing. 

The steady Dictions that are treated herein are assumed to give small 
deviations from the undisturbed flight path and the accelerated Dioticns 
are assumed to have small accelerations. Theoretical results based an 
this assumption for steady Dictions have, in general, been found to be 
reliable; however, the reliability of such results for unsteady motions 
is as yet unverified. The results presented herein cover a range of 
Mach number and aapect ratio greater than that for which the Mach line 
from the leading edge of the tip section intersects the trailing edge of 
the opposite tip section. 


SYMBOLS 


x,y rectangular coordinates (see fig. l) 

u 0 ,v Q induced flow velocities along x— and y-axes, respectively 

l coordinate in flight direction if this direction is inclined 

to x-axis 


u,v,w 


u 

w 


incremental flight velocities along x— , y— , and z-axes, 
respectively (see fig. 2) 

derivative of u with respect to time 

accelerated vertical Diotion 


Y 


undisturbed flight velocity 
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V* 


local flight Telocity after disturbance; used, to indicate 
inclination of flight direction to x— axis (see fig., l) 


p,q,r angular velocities about x— , y— , and. z— axes, respectively 
(see ‘fig. 2 ) 

a speed, of sound 

M stream Mach number (v/a) 


B 


= /m 2 - 


a 

a* 

a, 

t 

P 

c 

h 

b 

S 

Sw 

A 

Xcg 

p 

0 


Mach angle 


1 h) 


wing angle of attack in steady flight, .radians (w/7) 

local inclination of airfoil surface with respect to free 

stream, radians f =~ — ) 

3 \V + uy 

derivative of a with respect to time 
time following disturbance, seconds 
angle of sideslip, radians (v/v) 
chord 

wing semispan 
wing span 
total wing area 

region of integration over portion of wing surface (see fig. 3 ) 
aspect ratio 


(?-s) 


distance of origin of stability axes from the midchord point, 
measured along x-axis, positive ahead of midchord point 

masB density of air 

disturbance—' velocity potential on upper surface of airfoil 

auxiliary variables which replace x and y, respectively 
(see fig. l) 



k 


NACA TN No. 1706 


y a , q a indicates a transformation of origin of x— and y-axes or I— 
and Tj-axes from, leading edge of center section to leading 
edge of tip section (y a = y — & on right half -wing; 
y a = — y — h on left half-wing) 


AP 

ACp 


pressure difference "between lower and upper surfaces of air- 
foil, positive in direction of lift 

nondimensional coefficient expressing ratio of pressure 
difference "between lower and upper surfaces of airfoil 

to free— stream dynamic pressure 



X,Y,Z 


constant given "by equation ( 9 ) 

induced suction force cn wing tip per unit length of tip 

forces parallel to x— , y— , and z-axes, respectively (see 
fig. 2) 

longitudinal— force coefficient 





lateral— farce coefficient 



C Z 


vertical— force coefficient 




skin— friction drag coefficient 


rSkin-i friction drag' 


2 


L,M,N, moments about x— , y— , and z-axes, respectively (see fig. 2); 
M is also used to refer to Mach number 


rolling-moment coefficient 



( -m 


pitching-moment coefficient 
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°n 


yawing-moment coefficient 



Subscript: 


1,2 


contributions of normal pressures and skin friction, respectively. 


to 


also used to indicate component parts of 02 ^, Cy^ 


°%> °v 3111 Si 


Superscript : 


w,p contributions caused by vertical motion and rolling motion, 

respectively 


Whenever p, <p, r, 0, a, u, a,, and, u are used as subscripts, a 
nondimensional derivative is indicated and this derivative is the slope 


through zero. 


For example. 







i 


0 



TJhprlmed stability derivatives refer to principal body axes; 
primed stability derivatives refer to stability axes. 
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ANALTSIS 
General Concepts 


The coordinate axes and the symbols used in the' analysis of the 
rectangular wing are shown in figure 1. The derivation of the formulas 
for the surface— velocity-potential functions, pressure distributions, 
and stability derivatives is made initially with reference to principal 
body axes -which are fixed in the wing with the origin at the midchord of 

the center section ^-,0,0^. This system of axes is shown in figure 2(a). 

The transformation of these stability derivatives to a system of stability 
axes (fig. 2(b)) is discussed in the section entitled "Eesults and 
Discussion." 


The stability derivatives are determined from integrations of the 
forces and moments over the wing. .For vertical and pitching motions 
which yield equal and opposite suction forces along the edge of each 
wing tip, the only resultant forces and moments acting on the wing, if 
skin friction is neglected, are those caused by the pressures on the air- 
foil surfaces. These pressures are obtained from the familiar Bernoulli 
equation. 3n rolling, yawing, or sideslipping, however, unbalanced 
suction forces which produce lateral forces and yawing moments are 
induced along the wing tips in addition to the forces and moments 
resulting from the pressure normal to the wing surface. The subsequent 
analysis for the calculation of the stability derivatives is then 
resolved to a determination of the pressure distribution normal to the 
surface and the unbalanced suction forces along the wing-tip edges. 


The pressure 
(positive upward) 
general Bernoulli 


difference between the upper an d lower surfaces 
at any point on the wing is determined from the 
equation for small disturbances as 


AP = 2p 



30 30 'N 

3 1 3t / 


( 1 ) 


where Y* Is the local flight velocity and Z refers to a coordinate 
measured in the flight direction. The -term 30/3t expresses the effect 
of any unsteadiness in the flow. The velocity potential 0 in equa- 
tion (l) must be determined so as to satisfy the linearized partial- 
differential equation (with time dependency if the motion is unsteady) 
of the flow and the boundary conditions associated with the particular 
motion under consideration. Thus, the potential must give streamlines 
that are tangent to the airfoil surface and a pressure field that is 
continuous at all points exterior to the wing. Equation (l) shows that 
the pressure distribution on the wing is determined when the surface- 
potential function is found. 


The method of reference h is in general adaptable to the problem 
of obtaining the surface-potential function 0 in supersonic flight to 
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meet boundary conditions associated with, small steady motions, such ub 
vertical motions, rolling, sideslipping, and pitching. The method is 
an extension, to include tip effects, of the work of Puckett and others 
which uses the superposition of elementary source solutions to obtain 
the potential function. In cases where a point on the wing is influenced 
by two or more mutually interacting external fields, the interaction 
introduces difficulties In the solution for the surface potential. (See 
also reference 12.) If any point on the wing is influenced by only one 
independent external field, however, the potential function in a region 
affected by the wing tip may be obtained by integration of elementary 
source solutions solely over an appropriate area of the wing. The 
strength of these sources is shown to be a function only of the local 
slope of the airfoil surface with reference to the free— stream direction. 
Inasmuch as the slope of the airfoil surface with reference to the free— 
stream direction is known for a given motion, the distribution of sources 
is known and, consequently, the distribution of the surface-potential 
function is determined by an integration of the elementary source solu- 
tions over an appropriate area of the wing. 

As applied to the rectangular wing at supersonic speeds, the fore- 
going method of reference 4 for one independent external field is valid 
as long as the foremost Mach wave from one tip does not intersect the 
opposite tip, that is, for Mach numbers and aspect ratios for which 
AB 1. For this case, the potential at a point cn the top surface of 
a thin flat wing may be determined by means of equation (l4) of reference 4 
and is as follows: 


0(x,y a ) = 7. 


a 1 d£ dt] a 


I /(* - if - B 2 (y. - >i a ) 2 


( 2 ) 


where a* represents the local angle of at back of the airfoil surface 
at the point (1,%). Figure 3 shows a typical region Sy. for determining 
the potential at a point (x,y a ) in a rectangular wing. The figure 
shows the boundaries S„. over which the integration must be performed,- 
for a point (x,y a ) which is affected by the wing-tip region. If the 
point (x,y a ) is located at or inboard of the foremost Mach line from 
the tip, this point is unaffected by the tip region and is bounded 

by the leading edge and the Mach forecone from (x,y a ). Suppose that the 
surface potential 0(x,y) has been obtained from equation (2) or by 
.same other method, then the differentiation of 0 with respect to the 
coordinate in the free-stream direction determines the pressure distribu- 
tion by means of the Bernoulli relation, equation (l). 


Ths expressions for determining the surface potential and the pressure 
coefficient for unsteady motions are discussed in the section entitled 
"Derivation of Formulas . " 
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Derivation of Formulas 

The subsequent derivation of formulas for the various motions will, 
involve first the determination of distributions of surface potential and 
then the determination of surface— pressure distributions and any- 
unbalanced suction forces along the wing tips. The integrals required 
for these derivations and also those required for Hie stability derivatives 
are integrable either directly or after reduction by parts by means of 
standard formulas such as are given in reference 13 ; hence, the details 
for the integrations are not shown. In. the operations involving factoring 
from, radicals, care must be used to preserve the correct sign of the 
factors; for example, if 


3 r a < 0 


then 



For brevity, the final formulas are ammitted from the derivations and 
appear only in tables at the end of the paper. ThuB, the distributions 
of 0 and Acp are summarized in table I, and the stability derivatives 
are summarized in table H. 

* 

All the derivations are made specifically for a wing for which 
AB > 2, that is, for which the foremost Mach wave from a tip does not 
intersect the remote half— wing. The formulas in table I for the 
potential 0 and pressure coefficient Acp that are obtained for 

AB ^ 2 can be applied to wings in which 1 ^ AB < 2 by using the 
principle of symmetry and superposing separately each tip effect at the 
point under consideration to the value obtained for the infinitely long 
wing. A consideration of this superposition principle for the rectangular 
wing shows, however, that the stability derivatives which are obtained 
for AB ^ 2 apply as well to wings for which AB ^ 1 . A more detailed 
description of table II is given in the section entitled "Penults and 
Discussion. " 


Vertical, Pitching, and Longitudinal Motions 
Derivatives 8X1 ~^z ^ or pitching motion 

lateral axis through the midchord point, the local slope of the 
surface with respect to the free— stream direction is 


about a 
airfoil 


a 


t 


a + 



q. 


where a is the angle of attack in the absence of pitching. In order to 
obtain the potential distribution, this value of a* is substituted into 



NACA TN No. 1706 


9 


equation ( 2 ) and the double Integration for the "variables | and ^ 
is performed between the lim i ts indicated in figure 3 . The pressure 
coefficient is then obtained from equation (l) for steady motions as 


Ar.-n — 3. u 

P y o 


4 50 


( 3 ) 


These pressure coefficients are then differentiated with respect to a 
and q. The integrations of the respective distributions of Acp over 
the wing and conversion to nondimensional units then give the deriva- 
tives — and -C z . 

Derivatives CL and CL .-Derivatives and 0 are obtained 

~q “a ™q 

directly from, the pressure-coefficient distributions for angle of attack 
and pitching, thus 




where Ac p for angle of attack and pitching is given in table I. 

Derivatives C-^ and C^- .— At supersonic speeds the resultant 

pressure force on a rectangular wing of zero thickness acts normal to the 
surface as there is no suction at the leading edge. Thus, the forces 
in the x-direction arise solely from skin friction. On the assumption 

that the skin friction is Independent of a end q, the derivatives CV 

A a 

and C-y are zero. 


Derivatives 


-Cy C^, and — C-^.— The derivative -C^ is obtained 


from, the equation 


Then 


-Z = -a*C 2 ^ | (T + u) 2 S 


-C, 


Zu 


■ ▼ k h* c ^ (T + u)2 L. 
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The function — °frfc a i ne( 3- from table H. Its dependence on the 
incremental flight velocity TJ is indicated by giving B in the 

/ fy ^ ^ j2 -y 

form 1/ 1 and a* in the form . Then 

» a 2 Y + vl 


-0, -iL 

^ TBu 


W(v 


+ u 


> (- 


r(T + u) 2 


- 1 


V 


IcfjYW _ 


1 


u — >0 


The derivative 


Si 


is obtained from the equation 


Then 


M = | (V + u) 2 Sc 



The function C is obtained from table H, where its dependence 
•x 

on u is indicated by writing a’ and B in the same form given 
previously for the derivative — C^. Then 




% Y 


3A 


w(T + u) 

Rv + u ) 2 


.1 


- 1 


’u — ^0 


The derivative — < results from skin friction and is obtained 
from the equation 


— X = Cj^ §■ (Y + u) 2 S cos a 
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Than 



1 

Y Su 



(V + u) £ 


U- 


0 


Accelerated. Motions 

For accelerated motions in the vertical plane of symmetry, the 
pressure coefficient from equation (l) is 




GO 


The surface potential 0 in equation (4) for unsteady motions of thin 
airfoils in two-dimensional supersonic flow has "been derived in refer- 
ence l4. Jjl reference 10, the methods of references 4 and 12 for steady 
flow at supersonic speeds are extended in order to determine solutions 
for the surface potential and pressure coefficient for unsteady motions 
in three-dimensional flow. In the present analysis the solutions obtained 
in reference 10 will he utilized to calculate the derivatives in vertical 
motions with small constant accelerations. 


Derivatives -Cg. , C , and Cj^.— The surface potential 0 for 

uniformly accelerated motion as obtained for the region within the tip 
Mach cane is (reference 10, equation ( 31 ) ) 



(5) 


3h converting from the notation of reference 10 for a rectangular wing 
to the present notation, the following transformations are made: IJ = V, 

m = a, a = 0, 0 = B , c = J, ^ = 1, Ug = v w = — (x + y a B), u^ = — (x - y a B), 
du 2 Ac p 

dv w * P 2 


2B 
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la order to obtain 0 in idle region "between the tip Mach cones, y a in 
equation ( 5 ) is set equal to — J and, therefore. 


'-*r-SS5 (6) 

The pressure coefficient Acp contributed by the vertical accelerating 

motion is obtained by partial differentiation of 0 (equations ( 5 ) and ( 6 )) 
with respect to x and t, by letting t = 0 , and then by substituting 
these expressions for 80/ox and S0/dt in equation (4). This process 
yields in the region within the tip Mach ccnes 


A3 p (x,y a ) = 


ka x -1 /^aB 

5- COS l 

kB 2 ? B ' x 


+ 1) + 2B" 



and in the region between the tip Mach cones 

hrrrr . 

= (8) 

Equations (7) and ( 8 ) correspond to equation ( 33 ) of reference 10, 
after the appropriate transformations noted previously for 0 are made. 

The derivatives an d are then obtained by integration of 

the corresponding Acp-distributions and conversion to nondimensional 
units. The derivative Cj^ is shown to be zero by the use of assumptions 

similar to those noted previously for Cy^. 

Derivatives —I Cg^, C-m^, and — Cy^.— For small accelerations along the 

flight path, the potential will remain substantially unchanged. The 
increments in pressure caused by these accelerations, therefore, afe 

assumed to be negligible, and the derivatives -Cg^, Cj^, and are 

approximately zero . 


Bolling 


Derivative C-^.— In steady rolling motion with angular velocity p, 
the local slope of the airfoil surface with respect to the flow direction is 


a * = m _ 

a - T - 


P(% + B) 
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la order to obtain the potential distribution this value of a* is 
substituted into equation ( 2 ) and the double integration for the 
variables | and t] is performed between the limits indicated in 
figure 3 * The pressure coefficient is then obtained from equation ( 3 ). 
The derivative is obtained by integrating the mome nts of the 

Bernoulli pressure distribution for rolling given in table I and by 
converting this result to coefficient form. 


Derivatives 




and 



In a rolling motion. 


the lateral force 


and yawing moment relative to body axes result entirely from suction 
along the tips. These suction forces may be evaluated by applying a 
method suggested in reference 15 for incompressible flow and modified for 
compressibility effects in reference 3 . Thus, if the induced surface 
velocity normal to the wing tip is expressed as 


v 


o 


- + 


G 

i/-*k 


(9) 


where G is a constant, then the suction force per unit length of tip is 


Fg = otpG^ 

(A more general expression for edge suction that is still valid when the 
edge is inclined to the stream is given in reference 3 and recast in 
reference 7 - ) 

Consider the induced surface velocity normal to the tip of a wing 
rolling with an initial angle of attack a. This velocity is 


.v+p 30 /30V /3(A P 

T ° 3y \3y/ + \3y/ 

where the superscripts w and p refer to the potentials obtained for 
a vertical motion and a rolling motion, respectively. From table X 



^2y a B 



x 
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and partial differentiation of 


with respect to y yields 



( 10 ) 


where y a < 0. Very near the tip, — y a — ? 0 and 



( 11 ) 


The potential in rolling 0^ is given in table I. By partial differentia- 
tion of 0P with respect to y and then by letting — y a — *0, there 
results 




( 12 )' 


The resultant induced surface velocity normal to the wing tip as 
— y„ — =>0 is obtained by adding equations (ll) and {12.). Thus 

CL 



2 fx 

jt ys 



± Ya 


(13) 


The plus signs before p and Y refer to the right wing tip and the 
negative signs refer to the left wing tip. 

Yery near the wing tip, equation (13) has the same form as equa- 
tion ( 9 ) and, therefore, the total suction force per unit length along 
the wing tip is 


F 


S = 


jtpG^ = 





(14) 


In equation (l4) only the term ±- 


8pxpYa 
itB ■ 



will 


give rise 


to a lateral force and a yawing moment which 
this term along the wing tips. These forces 
converted to nondimensional form to give the 


are obtained by integrating 

and moments are then 

derivatives CV and CL 
X P Op 
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Sideslip 

The pressure coefficient obtained from equation (l) for steady- 
flight is 



where V* and l are measured in the flight direction. If sideslip 
occurs the flight direction is inclined relative to the x-axis by the 
sideslip angle p. The rectangular wing in sideslip, therefore, becomes 
equivalent to a yawed wing with the leading wing tip raked out and 
the trailing wing tip raked in. If the Kutta-Joukowski condition at the 
trailing wing tip is neglected, the potential function for the yawed 
rectangular plan form may be obtained by the method of reference 4. In 
reference 11, the method of reference If- is extended in order to obtain 
solutions for edges for which the Kutta-Joukowski requirement must be 
satisfied. 

Physical considerations suggest, however, that for small sideslip 
angles, the actual flow for typical ! y rounded wing tips would in general 
be unlikely to conform to the Kutta-Joukcwski condition along the trailing 
wing tip. The edge suction for a lifting wing arises because of the flow 
from the bottom surface to the top surface around the side edge. This 
flow may be presumed to go around any boundary layer that may be present. 
The local boundary layer thus experiences the edge suction. Rough 
calculations suggest that the edge suction per unit area is approximately 
constant from the leading edge to the point of maximum profile thickness, 
and then increases rapidly from the- point of maximum thickness to the 
trailing edge. The pressure gradient is therefore favorable and the 
flow at the side edge is pot expected to separate. This condition should 
persist for smsl 1 or moderate amounts of sideslip until the additional 
pressure increment caused by sideslip produces a strong' adverse pressure 
gradient. Further theoretical and experimental investigation is required 
to obtain quantitative results regard i ng these phenomena. On the basis of 
the foregoing considerations, it will be assumed in the present analysis 
that the Kutta-Joukcwski condition is not satisfied along the trailing 
wing tip. The effect of satisfying the Kutta-Joukowski condition along 
the trailing wing tip in sideslip is discussed in this analysis and also 
in the section entitled "Results and Discussion. " 


Derivative 



— The potential corresponding to a thin rectangular 


wing at an angle of attack and a finite angle of sideslip may be obtained 
from reference 4, equation (20). The corresponding pressure distribution 
may be obtained from reference 11, appendix C, equation (c4). These 
solutions from references 4 and 11 were simplified to the approximate 
form for sma.1 1 angles of sideslip (p^ « l) and converted to the present 
notation with respect to axes shown in figure 1. The distributions 
for 0 and Acp caused by combined vertical motion and sideslip are 
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given in table I. The regions for -which these expressions for p a n d Acp 

are applicable are bounded by Mach lines with respect to the stream, 
velocity V* -which is inclined to the z-arl b by the sideslip angle p. 

As noted previously, these expressions do not satisfy the Kirtrba-Jaukawski 
condition along the trailing wing tip. As indicated in reference 11, 
however, the Kutta-J oukovskL condition along the trailing wing tip merely 
cancels the radical term in the expression for dcp within the Mach cane 
from the trailing wing tip. 

A consideration of the foregoing Acp-distributians Indicates that 

as a result of sideslip the lift within the Mach cone from the leading 
wing tip is increased, whereas the lift within the Mach cone from the 
trailing wing tip is decreased. A rolling moment is thereby produced. 
Furthermore, as a result of sideslip, the Mach lines are shifted toward 
the trailing wing tip, and this shift contributes an additional ro l ling 
moment ; The magnitude of the rolling moment caused by sideslip iB given 
in table H in terms of the ncnd imenslanal derivative /C 7 A 

V ^0—0 


Derivatives C-, 


f and C_ .— The derivatives 

3 ~fL 


°Ti 


P 


and C. 


*P 


can result 


solely from, suction forces which are induced at the wing tips. These suction 
forces for sideslipping motion were evaluated by a method similar to that 
described previously for obtain 0^ and V Tbs treats for sideslip 

was based on the conclusion, noted previously, that the Kut t a-J oukowski 
condition is unlikely to be satisfied for typically rounded wing tips at 
smal l angles of sideslip. The potential p for determining the induced 
velocity normal to the wing tip was obtained from table I. The resultant 
lateral force and yawing moment are given in nondimenslonal form in table II. 


Tawing 

In yawing flight, the stream velocity varies linearly along the span. 
This effect introduces variations of both dynamic pressure and compressi- 
bility effects along the wing span. The surface potential as expressed in 
equation ( 2 ) satisfies the linearized potential equation for a uniform 
stream Mach number, but is inadequate to account for the compressibil i ty 
effects associated with a spanwise variation of stream Mach number. (See 
reference 7. ) The caBe of the trapezoidal, wing with tips cut off alon g 
the Mach lines (raked tips) was analyzed in reference 7* It was shown 
that the pressure distribution could be obtained by application of the 
Ackeret two-dimensional theory modified by using the local Mach num ber at 
each spanwise station as affected by the yawing. Inclusion of the span- 
wise variation in Mach number was demonstrated to have a profound effect 
on the pressure distribution. 
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The addition of suitable triangular tips to the aforementioned 
trapezoidal wing converts it into a rectangular wing. The added tips 
lie wholly within the tip Mach cones and thus their addition does not 
alter the pressures on the trapezoidal .portions. A rigorous solution 
for the pressures on the tip portions cannot yet be demonstrated. However, 
an expression that appears plausible has been obtained. This pressure 
distribution for the tip portions is derived by superposing on the 
Ackeret pressure distribution, as modified by local Mach number, an 
appropriate' function which fulfills the boundary condition for no 
pressure discontinuities in the region exterior to the wing. This function 
thus represents the effect of the wing cut-off and is designated herein 
as the tip effect. The pressure difference AP at any point according 
to the Ackeret theory based an local Mach number is (reference 7 ): 



Equation ( 15 ) shows that -the pressure distribution for an Infinitely 
long wing which has a steady yawing velocity r and vertical velocity w 
is expressed by two components. One of these components is proportional 
to w, is constant, and gives the pressure distribution contributed by 
an angle of attack in straight flight. The other component is proporti on al 
to wr, gives a linear antlsymmetrical distribution with respect to y, 
and expresses the pressure distribution contributed by yawing. 

o 

It will be recalled that the solution for steady rolling, treated in 
a preceding section, resulted likewise in a pressure distribution 
proportional to y in the region between the tip Mach cones. The pressure 
distributions contributed by rolling and by yawing are thus proportional 
in the region between the tip Mach cones. The wing cut-off is effected by 
canceling the disturbance pressures outboard of the desired tip location 
by means of a function that satisfies the boundary conditions on the wing. 
Because the two pressure distributions to be canceled correspond in the 
yawing and rolling cases, the incremental pressure function or tip effect 
for each case evidently must reduce to forms which will have the same 
factor of proportionality in the entire plane of the wing outboard of 
the tip. It seems reasonable to assume, therefore, that for small yawing 
motions the two pressure distributions will also have very nearly the 
same factor of proportionality within the tip Mach cones. 

The proportionality constant between the pressure distributions for 
rolling and yawing motions may be determined by a comparison of the cases 
of rolling and, yawing in column 4 of table I. The pressure coefficient 
per unit yawing velocity is seen to be a/B^ times the pressure coeffi- 
cient per unit rolling velocity, or 
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(Ac) =^-fAcp) (16) 

V ■^'Yawing B 2 p V. V Rolling 

where equation (l6) will apply over the whole wing. 

Derivative C .— The preceding analysis indicated' that the pressure 
_ hi 

» distribution per unit yawing velocity is in a simple ratio to that 
produced per unit ro ll ing velocity (equation (l6)). Thus 



The derivative C 7 has "been derived previously and is given in table H. 
P 


Derivatives 


CL. and G .— When the wing yaws, the antiBymmetrical 
x r wr 


pressure distribution which is indicated by equation (15) will produce 
unbalanced suction forces at the right and left wing tips and thereby 
give rise to lateral forces and yawing moments. In addition, skin friction 
will contribute a yawing moment. 


It appears that a reasonable although approximate evaluation of the 
tip suction forces in yawing can be obtained by means of the correspondence 
of yawing with rolling as utilized previously in deriving equation (l6). 
This procedure does not satisfy the Kutta-Joukowski requirement in the 
sideslip component of the stream velocity in yawing; however, this 
theoretical deviation is likely to be very small in the actual flow. On 
the basis of these considerations, the induced suction forces an the 
wing tips per unit yawing velocity will be related in the ratio — to 

B 2 

the corresponding induced suction forces per unit rolling velocity which 


were derived previously ^section entitled "Derivatives Cy^ and C^'). 
The contributions of the tip suction forces to Cy and C n ^ are, 
therefore. 


and 





a 

B 2 



where Cy and C 

P P 


are given in table H. 
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The effect of skin friction on the yawing moment due to yaw ing is 


Ng = cos a 



J -h. u 0 ^0 2 


(V 


- ry) 2 + r Q - x) (x - 3 + y 


dx dy 


where the first bracketed term expresses the square of the resultant 
local velocity and f3 is the local angle of sideslip: 


y 

Eliminating second— order terms and terms corresponding to symmetrical 
drag forces and converting Eg to coefficient form yields 



%_ 

4ch^ 




dx dy 


RESULTS AND DISCUSSION 


As noted in the preceding analysis, the n nnd l-merm l mml stability 
derivatives which are presented in table H were derived with reference 

to principal body axes with the origin at point (§,°,o). These results 

may be transformed by means of the equations in the last column of 
table H to apply to stability axes with the origin at an arbitrary 
distance x Q g from the midchord point. The stability axes are shown 

in figure 2(b) and are obtained by a rotation of the principal body 
axes (fig. 2(a)) through an angle a; the origin is then shifted a 
distance %ca along the new x-axis. The conversion to stability axes 
was obtained by means of the transformation fo rmal a s given in reference 16 , 
with the omission of relatively unimp ortant terms c omp ared to uni ty, 
such as a 2 . 

The formulas for the derivatives given in table H with reference to 

principal body axes are shown plotted in figures 4 and 5 against the 

( ' 0 

Derivative — C-^ and those derivatives equal to zero 

are omitted from the. figures^ These curves. show the variation of the 

stability derivatives with aspect ratio for constant Mach number. The 
variation with Mach number for constant aspect ratio is not directly 
indicated, although it can be determined from the curves. These data 
are shown in figure 4 for the lateral stability derivatives and in 
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figure 5 for the longitudinal stability derivatives. The data in 
figures 4 and 5 are to he used in conjunction -with the transformation 
formulas presented in table II .to evaluate the derivatives with respect 
to stability axes. In the evaluation of these derivatives, many of the 
terms are likely to be relatively small; therefore, the expressions will be 
noticeably simplified -when such terms are neglected in the computations. 

' The results of the present investigation have 'been derived on the 
assumptions of zero thickness and small disturbances. Potential flow is 
assumed except in the case of C^, and in which sk in friction is 

considered. The practical effects of the Kutta— Joukowski requirements 
■which are introduced at the wing tips in sideslip and yawing are not 
definitely known. On this basis, the data shown in figures 4 and 5 are 
expected to apply in getneral to thin wing sections for small steady 
motions, motions with small accelerations, or oscillatory motions of 
low frequency in which boundary— layer effects are not expected to be 
important. The applicability of the present theory to Mach numbers in 
the vicinity of unity, very high Mach numbers, or for very low aspect 
ratios is -uncertain. 

The data in figure 4(a) show that at supersonic speeds the 
b3c z 

derivative — is negative in contrast with the behavior at subsonic 

a 

speeds where positive values would be obtained. This phenomenon was 
pointed out for the infinitely long wing in reference 7 and its physical 
significance elaborated upon. For stability axes, the formula for C 2 r * 

(table H) indicates that another reversal in sign to a positive value 
occurs as the Mach number is increased beyond approximately 1.4l for 
typical rectangular wings. (Also see reference 7 for the infinitely 
long wing. ) 


The suction force at the leading edge of rectangular wings vanishes 
at supersonic speeds. This factor should have an important influence an 
the derivatives Qq and as supersonic speeds are attained. In 


S' 


p Pc 

(stability axes), the results of the present analysis 


the case of 

indicate that at supersonic speeds the sign of 


S' 


will have positive 


values in many typical cases in contrast to negative values normally 
obtained at subsonic speeds. In the case of or Cq^.*, the loss of 

leading-edge suction tends to be compensated by the spanwise compressi- 
bility effects associated with supersonic speeds. 


As noted previously in the analysis, the Eitta— Joukowski condition 
is unlikely to be satisfied along the trailing wing tip for a typically 
rounded wing tip at small angles of sideslip. Therefore, the results 
for Ci |n figure 4(a) are applicable where the Eutta-J oukovski 
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condition along the wing trailing edge is not satisfied. In order to 
determine the effect on C-^ of satisfying the Eu.t t a-J oukovski condition 

along the trailing wing tip, the formula for C 2 ^ 'which meets this 


requirement was also obtained and, is as follows: 



aA _ 3 + 2 B 2 \ 

3A 2 B 2 / 


A comparison of this formula with the data for 


C-j given in figure 4(a) 
P 


indicates that the effect of satisfying the Eutta-Joukcwski condition 
along the trailing wing tip reduces negatively the values of from 

p 

those obtained by neglecting the Eut ta-J oukowski condition. For example 
for B = 1 and AB = 4, -when the Eutta-Joukcwski condition along the 
trailing wing tip is neglected, C|p = -0 . 083 a; and -when the Eutta-Joukcwski 

condition is satisfied along the trailing w ing tip, = 0.l46a. Thus, 

it is expected that when the sideslip angle becomes large, the dihedral 
effect should be reduced significantly because of the Eutta— 

Joukowski condition along idle trailing wing tip. 


The longitudinal, stability derivatives in figure 5 refer to an 
axis whose origin is located at the midchord point. The data in figure 5(b) 
for BCm show that rect ang ular wings, with reference to this origin, 
have an increasingly unstable pitching moment with decreasing aspect 
ratio which corresponds to a forward shift in the aerodynamic center. 

For infinite aspect ratio, the aerodynamic center is located at the 
midchord point or BCL =0. If the aspect ratio is decreased to a 


value of 4 for a Mach number of 1.4l, figure 5 indicates a forward 
shift of the aerodynamic center of 0.025 chord. With constant Mach 


number, the ratio B< 




is obtained from figure 5 solely as a 


function of AB. These data indicate that with constant aspect ratio 
and increasing Mach number, the aerodynamic center will shift rearward. 
For an aspect ratio of 4, an increase in Mach number from 1.4 to 1.9 will 
shift the aerodynamic center rearward 0.01 chord. 


The derivative — given in table H for inf ini te aspect ratio 

is negative which indicates negative damping or instability. Tha 

ratio -Cjjj^y-C ^ from table H gives the location of the center of pressure 

of the resultant lift contributed by <x. By taking this ratio for 
infinite aspect ration the center of pressure is found to be located at 

a point c behind the leading edge. The negative dainplng produced 
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"by a, therefore , gives 
gravity locations ahead 
"by a are minimi zed by 

due to Cjr. disappears 


an unstable pitching moment for center-of— 
of x = j c. These unstable tendencies caused 
the effects of finite span and the instability 
s M +1 

entirely if A Bn . 

3 


CfflJCLOSIQNS 


A theoretical investigation has been made by means of the linearized 
theory to obtain formulas for the surface-velocity-potential functions ■, 
surface— pressure distributions , and, stability derivatives for various 
motions at supersonic speeds for rectangular wings of zero thickness 
without dihedral. The investigation included steady and accelerating 
vertical and longitudinal motions and steady rolling, yawing, sideslipping, 
and pitching for Mach numbers and aspect ratios greater than those for 
which the Mach line from the leading edge of the tip section intersects 
the trailing edge of the opposite tip section. 

The following significant conclusions have been obtained for this 
investigation: 

1. At supersonic speeds for Mach numbers smaller than approximately 
1.4l, positive yawing generally results in a negative rolling moment in 
contrast to the behavior at subsonic speeds where a positive rolling 
moment is produced. 

2. The attainment of supersonic speed produces a significant change 
in the positive direction of the yawing moment per unit rolling velocity. 


3 . For Infinite aspect ratio, a constant vertical acceleration 
causes a negative damping in the vertical motion, and, an unstable pitching 

p 

moment for center— of-gravity locations ahead of the ^ —chord point. 

These unstable tendencies are mini ml. zed by the effects of finite span 
and the instability due to the rate of change of lift with vertical accelera- 
tion disappears entirely if A\ju 2 - 1 ^ M 2 - + - 1 where A is the aspect 
ratio and M is the Mach number. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., June 30* 19^8 
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Wing tip 


V 



Figure 3.- Region of integration for obtaining supersonic velocity -potential 
function for rectangular wing of finite wing span. Region of integration 
for point (x,y a ) shown brosshatched; y a = y -’h on right half -wing; 

y a = -y - h on left half -wing. 
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(a) Vertical -force -coefficient 
derivatives. 


(b) Pitching -moment -coefficient 
derivatives. 


Figure 5.- Variation of Supersonic longitudinal stability derivatives with aspect 
ratio -Mach number parameter. Derivatives with respect to principal body 
axes; thin flat rectangular wing. (See table H for conversion to stability 
axes.) 
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